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The Euclidean action of non-Abelian gauge theories with adjoint dynamical charges 
(gluons or gluinos) at non-zero temperature T is invariant against topologically 
non-trivial gauge transformations in the Z(JV) C center of the SU(N) gauge group. 
The Polyakov loop measures the free energy of fundamental static charges (in- 
finitely heavy test quarks) and is an order parameter for the spontaneous break- 
down of the center symmetry. In SU(N) Yang-Mills theory the Z(7V) C symmetry 
is unbroken in the low-temperature confined phase and spontaneously broken in 
the high-temperature deconfincd phase. In 4-dimensional SU(2) Yang-Mills the- 
ory the deconfinement phase transition is of second order and is in the universality 
class of the 3-dimensional Ising model. In the SU(3) theory, on the other hand, 
the transition is first order and its bulk physics is not universal. When a chemical 
potential fi is used to generate a non-zero baryon density of test quarks, the first 
order deconfinement transition line extends into the (fi, T)-planc. It terminates at 
a critical endpoint which also is in the universality class of the 3-dimensional Ising 
model. At a first order phase transition the confined and deconfincd phases coexist 
and are separated by confined-deconfined interfaces. Similarly, the three distinct 
high-temperature phases of SU (3) Yang-Mills theory are separated by deconfined- 
deconfined domain walls. As one approaches the deconfinement phase transition 
from the high-temperature side, a deconfined-deconfined domain wall splits into a 
pair of confined-deconfined interfaces and becomes completely wet by the confined 
phase. Complete wetting is a universal interface phenomenon that arises despite 
the fact that the bulk physics is non-universal. In supersymmetric SU(3) Yang- 
Mills theory, a Z(3) x chiral symmetry is spontaneously broken in the confined 
phase and restored in the deconfincd phase. As one approaches the deconfinement 
phase transition from the low-temperature side, a confined-confined domain wall 
splits into a pair of confined-deconfined interfaces and thus becomes completely 
wet by the deconfined phase. This allows a confining string to end on a confined- 
confined domain wall as first suggested by Witten based on M-theory. Deconfincd 
gluons and static test quarks are sensitive to spatial boundary conditions. For 
example, on a periodic torus the Gauss law forbids the existence of a single static 
quark. On the other hand, on a C-periodic torus (which is periodic up to charge 
conjugation) a single static quark can exist. As a paradoxical consequence of the 
presence of deconfined-deconfined domain walls, in very long C-periodic cylinders 
quarks are "confined" even in the deconfined phase. 
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1 Introduction and Summary 

Understanding the non-perturbative dynamics of quarks and gluons is a highly 
non-trivial task. Since at low temperatures these particles are confined into 
hadrons, the quark and gluon fields that appear in the QCD Lagrangian do 
not provide direct insight into the non-perturbative physics of confinement. In 
the real world with light quarks confinement is accompanied by spontaneous 
chiral symmetry breaking and the presence of light pseudo-Goldstone pions. 
The pion fields that appear in the effective chiral Lagrangian do provide direct 
physical insight into the low-energy physics of QCD, but they do not shed light 
on the dynamical mechanism that is responsible for chiral symmetry breaking 
or confinement itself. 

In order to gain some insight into the physics of confinement, it is very 
interesting to study QCD in extreme conditions of high temperature and large 
baryon chemical potential. Due to asymptotic freedom, confinement is then 
lost and we can hope to learn more about confinement by understanding how 
it turns into deconfinement. For the theory with light quarks this question is 
still very difficult to address from first principles, although a lot of interesting 
results have been obtained using lattice field theory as well as phenomenological 
models. 

Here we decide to study a hypothetical world without dynamical quarks. 
Then the only dynamical color charge carriers are the gluons. Quarks appear 
only as infinitely heavy test charges that allow us to probe the gluon dynamics. 
In the resulting Yang-Mills theory confined and deconfined phases are easier 
to understand than in full QCD because they are characterized by distinct 
symmetry properties. As first pointed out by 't Hooft El, the symmetry that 
distinguishes confinement from deconfinement is the Z(7V) C center symmetry 
of the SU(N) gauge group. 

The selection of topics covered in this article is strongly biased by the 
authors' personal interests. Still, we hope that it will give a useful overview 
of the many interesting phenomena related to the center symmetry and its 
spontaneous breakdown. There is an extensive literature on the subject and 
we will not be able to discuss all important contributions. We like to apologize 
to those whose work will not be mentioned here. 

1.1 Static Quarks and Dynamical Gluons 

Gluons transform in the adjoint representation of SU(N) and are Z(iV) c neu- 
tral. A single test quark transforms in the fundamental representation and 
carries a unit of ~L(N) C charge, while N test quarks (forming a "test baryon") 
are together Z(7V) C neutral. In Yang-Mills theory a confined phase is simply 
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characterized by the fact that a non-zero Z(JV) C charge (for example, a single 
test quark) costs an infinite amount of free energy. In a deconfined phase, on 
the other hand, the free energy of a single static quark is finite. In full QCD 
this simple confinement criterion fails because in the presence of dynamical 
quarks the free energy of a single static test quark is always finite. This is 
because the color flux string emanating from the static quark can end in a dy- 
namical anti-quark. The static-quark-dynamical-anti-quark meson has a finite 
free energy. In Yang-Mills theory, on the other hand, there are no dynamical 
quarks or anti-quarks and the string emanating from a test quark is absolutely 
unbreakable. Since it carries the X(N) C flux of the test quark, the string can- 
not end in a dynamical gluon which is ~L{N) C neutral. Instead it must go to 
infinity. Since the confining string has a non-zero tension (free energy per unit 
length) the static quark with the infinite flux string emanating from it costs 
an infinite amount of free energy. 

Polyakovi and Susskindi were first to point out that the operator that 
describes a static quark is a Wilson loop closed around the periodic Euclidean 
time direction — the so-called Polyakov loop <E>. The free energy F of the test 
quark is determined by the thermal expectation value of the Polyakov loop 
as ($) = exp(— (3F). Hence, in a confined phase in which F — oo one has 
($) = 0. This can be understood as follows. At low temperatures T the extent 
(3 = 1/T of the Euclidean time direction is large. Then the Polyakov loop 
extends through many space-time regions that are essentially uncorrelated in 
color space. As a consequence, the Polyakov loop behaves essentially as a 
random variable and its thermal expectation value averages to zero. In the 
high-temperature deconfined phase, on the other hand, F is finite and hence 
(<I>) ^ 0. In that case the Euclidean time direction is short and the Polyakov 
loop extends only through a small space-time region that is highly correlated in 
color space. Then the Polyakov loop orders and picks up a non-zero expectation 
value. 

Under a twisted gauge transformation that is periodic up to an element z 
of the center 7.(N) C , the Polyakov loop changes by a factor z, while the action 
remains invariant. Consequently, the Polyakov loop is an order parameter for 
center symmetry breaking. A non-zero value of ($) implies that the center 
symmetry is spontaneously broken. This is the case in the high-temperature 
deconfined phase. In the confined phase, on the other hand, the center sym- 
metry is unbroken. This result is counter-intuitive, because one might expect 
spontaneous symmetry breaking to occur only at low temperatures. However, 
it should also not be too surprising because, for example, even the Ising model 
displays spontaneous symmetry breaking at high temperature when it is formu- 
lated in terms of dual variables. It should be noted that the center symmetry 
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is absent in full QCD with dynamical quarks because their presence leads to 
explicit Z(3) c breaking. 

1.2 Phase Transitions and Interface Dynamics 

As the order parameter that distinguishes confinement from dcconfincmcnt, the 
Polyakov loop plays a central role in understanding the physics of the high- 
temperature phase transition. In fact, just as the pion field of the effective 
chiral Lagrangian provides direct insight into the low-energy physics of QCD, 
an effective action for the Polyakov loop provides direct insight into the high- 
temperature physics close to the confinement-deconfinement phase transition. 
In the case of SU (2) Yang-Mills theory, the Polyakov loop is represented by a 
real-valued scalar field in three dimensions with a Z(2) c symmetry. Hence, the 
effective theory for the Polyakov loop is a 3-d $ 4 theory. This theory is in the 
universality class of the 3-d Ising model. It was first noted by McLerran and 
Svetitskyta and by Svetitsky and YaffeU that the second order phase transition 
in the 3-d Ising model has the same universal features as the deconfinement 
phase transition in the SU(2) Yang-Mills theory. The first numerical lattice 
simulations indicating a second order SU(2) deconfinement phase transition 
were performed inl3. A detailed analysis of the critical exponents shows that 
they are indeed the same as in the 3-d Ising model. The lattice calculations 
of the value of the critical temperature have been extrapolated reliably to the 
continuum limit 13. 

In SU (3) Yang-Mills theory the Polyakov loop is described by a complex- 
valued field $ in three dimensions. The effective theory is again a 3-d <E> 4 
theory but now with a Z(3) c symmetry. A lattice model with the same sym- 
metries is the 3-A. 3-state Potts model □ which is known to have a first order 
phase transition!!!]. In SU(3) Yang-Mills theory the_deconfinement phase tran- 
sition is therefore also expected to be first order so that the bulk physics 
is not universal. Lattice gauge theory has indeed provided detailecL-numeri- 
cal evidence for a first order transition in SU{3) Yang-Mills theory E3. At a 
first order phase transition both phases coexist with each other and are sep- 
arated by interfaces. For example, at the deconfinement phase transition the 
low-temperature confined phase and the high-temperature deconfined phase 
are separated by confined-deconfined interfaces. The simplest mechanism for 
the dynamics of a first order phase transition is bubble nucleation. For ex- 
ample, in boiling water bubbles of the gas phase form and grow in the liquid 
phase. Similarly, when one heats up a system of glueballs (confined gluons) 
beyond the temperature of the deconfinement phase transition, one expects 
the formation of bubbles of gluon plasma (deconfined gluons) which grow and 
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ultimately fill the entire volume with deconfined phase. Due to the Z(3) c 
symmetry, there are in fact three distinct deconfined phases. Hence, bubbles 
of different deconfined phases may nucleate independently in different parts 
of the volume. When these bubbles expand and collide they cannot coalesce 
and form a single bubble because they contain distinct types of deconfined 
phase. Consequently, when the phase transition is completed, we are left with 
different deconfined regions separated by deconfined-deconfined domain walls. 
These domain walls can exist at all temperatures above the phase transition 
point, while confined-deconfincd interfaces are stable only at the phase transi- 
tion. Deconfined-deconfined domain walls have a non-zero tension (free energy 
per unit area) . Due to asymptotic freedom the domain wall tension can-be cal- 
culated analytically at high temperatures using semi-classical methods EJ. The 
confined-deconfined, interface tension, on the other hand, is accessible only to 
lattice simulationsli-1 because the phase transition itself occurs at a temperature 
at which analytic methods based on asymptotic freedom are not applicable. It 
should be noted that in full QCD deconfined domains are unstable or at most 
meta-stable, because dynamical quarks break the center symmetry explicitly. 
The fate p£ a meta-stable deconfined domain in the early universe has been 
discussed Ej. Via a tunneling process it turns into the stable deconfined phase 
at temperatures around 10 TeV, i.e. much before it reaches the deconfinement 
phase transition at temperatures around 100 MeV. 

When one cools the gluon plasma down below the phase transition tem- 
perature, one expects bubbles of the confined phase to nucleate in the bulk 
deconfined phase. These bubbles grow and ultimately fill the entire volume 
with confined (glueball) phase. Interestingly, in the presence of deconfined- 
deconfined domain walls there is a competing mechanism for generating the 
confined phase. Rather than nucleating in the deconfined bulk (which costs 
a large amount of confined-deconfined interface tension), it is energetically fa- 
vorable for the confined phase to appear at a deconfined-deconfined domain 
siall, thus making use of the free energy provided by the domain wall tension 
E3. In fact, the confined phase appears as a thin film that completely covers 
the domain wall. Hence, the deconfined-deconfined domain wall splits into 
a pair of confined-deconfined interfaces. Such phenomena are well known in 
condensed matter physics under the name of complete wetting. For example, 
in the human eye a film of liquid of tears completely wets the eye and thus 
separates it from the surrounding air. The alternative to complete wetting is 
incomplete wetting. Then lens shaped droplets of one phase appear at the 
domain wall separating two other phases. Interestingly, complete wetting is a 
universal phenomenon that occurs despite the fact that the phase transition 
is first order and the bulk physics is non-universal. In fact, complete wetting 
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cannot occur at a second order bulk phase transition because it requires inter- 
faces and thus phase coexistence. Complete wetting is characterized by a set 
of universal critical exponents which can be extracted from the effective 3-d 
Z(3) c symmetric $ 4 theory for the Polyakov loopE3. 

The Z(3) c center symmetry also plays an important role in the high- 
temperature physics of supersymmetric Yang-Mills theory. The supersym- 
metry partners of the gluons are adjoint Weyl fermions — the gluinos. Like 
gluons, gluinos are Z(3) c neutral and their confinement or deconfinement can 
again be characterized by the value of the Polyakov loop. In particular, the 
Z(3) c flux of an external fundamental test quark cannot end at a dynamical 
adjoint gluino and must hence go to infinity. In the confined phase this still 
costs an infinite amount of free energy. Besides the Z(3) c symmetry, super- 
symmetric Yang-Mills theory has a Z(3) x chiral symmetry. This symmetry 
is a remnant of the U(1)r symmetry that is explicitly broken down to Z(3) x 
by an anomaly. In the low-temperature confined phase the Z(3) x symmetry 
is spontaneously broken by the formation of a gluino condensate, while the 
Z(3) c symmetry remains unbroken. At high temperatures one expects chiral 
symmetry to be restored and the Z(3) c symmetry to be spontaneously broken 
due to deconfinement. 

Due to the spontaneous breakdown of the Z(3) x chiral symmetry there 
are three distinct confined phases in supersymmetric SU(S) Yang-Mills theory 
that are distinguished by the value of the gluino condensate. Below the decon- 
finement phase transition these phases coexist and are separated by confined- 
confined domain walls. When such a domain wall is heated up to the phase 
transition temperature, one of the three high-temperature deconfined phases 
forms a complete wetting layer and the copfined-confined domain wall splits 
into a pair of confined-deconfined interfaces!^. The fact that the wetting layer 
consists of deconfined phase has interesting consequences. It allows a color flux 
string emanating from a static test quark to end on the wall because inside the 
deconfined wetting layer the flux spreads out and no longer costs infinite en- 
ergy when transported to infinity. The phenomenon of strings ending on walls 
in supersymmetric. Yang-Mills theory was first predicted by Witten in the con- 
text of M-theoryEl The color flux string then appears as a fundamental string 
and the confined-confined domain wall manifests itself as a D-brane — a nat- 
ural place for strings to end. Complete wetting provides a field theoretical 
explanation of this phenomenon. 
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1.3 The Center — A Subtle Symmetry 

It should be noted that — unlike other symmetries — the center symmetry is 
not a symmetry of the Hamiltonian which describes the time evolution of the 
system. Instead, it is a symmetry of the transfer matrix which describes the 
spatial evolution of the system in thermal equilibrium. This subtle difference 
is at the origin of a controversy about the physical reality of distinct high- 
temperature phases and of deconfined-deconfined domain walls. It has been 
argued that these domain walls are just Euclidean field ..configurations (like 
instantons) which play no role in Minkowski space-time c3. For example, in 
contrast to other topological excitations like monopoles or cosmic strings, one 
cannot expect a deconfined-deconfined domain wall to appear as an eigenstate 
of the theory's Hamiltonian. This should not be too surprising because the 
deconfincd phases themselves are thermal mixtures of pure states weighted by 
their Boltzmann factors. Of course, the entire spectrum of the Hamiltonian is 
completely independent of temperature. The spectrum of the transfer matrix, 
on the other hand, is very sensitive to the temperature which controls the size 
of the Euclidean time dimension. Despite the fact that deconfined-deconfined 
domain walls are not energy eigenstates, their presence in the thermally equi- 
librated gluon plasma has observable consequences — at least in Gedanken 
experiments.^ For example, the free energy of a static test quark in a large box 
filled with gluon plasma depends crucially on the shape of the box. This effect 
is due to the presence of distinct high-temperature phases and the deconfined- 
deconfined domain walls that separate them. 

To study the free energy of a single static quark in a finite box one must 
choose appropriate boundary conditions. Standard periodic boundary condi- 
tions are inappropriate because — as a consequence of the Z(3) c Gauss law — 
a single quark cannot even exist on a torus E3. In a periodic volume the color 
flux emanating from a quark cannot go to infinity and must thus end in an 
anti-quark. The system as a whole is thus Z(3) c neutral. In a finite periodic 
volume the expectation value of the Polyakov loop always vanishes even in the 
deconfined phase. As a consequence of the Z(3) c symmetry there is tunnel- 
ing between the three deconfined phases which averages the Polyakov loop to 
zero. Irrespective of the dynamics, a single quark in a periodic box always has 
infinite free energy just for topological reasons. An appropriate alternative are 
.Cnperiodic boundary conditions which are periodic up to charge conjugation 
r'lr 4 !. In that case, the C-periodic copy of a quark is an anti-quark outside the 
box that can absorb the flux. The C-periodic box itself contains just a single 

a In the real world the presence of dynamical quarks leads to an explicit Z(3) c breaking 
which eliminates two of the three deconfined phases as well as the deconfined-deconfined 
domain walls. 
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quark (or more precisely a superposition of a quark and an anti-quark) and is 
thus charged. C-periodic boundary conditions lead to an explicit breaking of 
the center symmetry that disappears only in the infinite volume limited. As a 
consequence, the Polyakov loop no longer averages to zero and the free energy 
of a static quark becomes finite. In the confined phase the free energy goes to 
infinity in the infinite volume limit, while it remains finite in the deconfined 
phase. C-periodic boundary conditions are useful in a Gedanken experiment 
that demonstrates the physically observable consequences of the presence of 
deconfined-deconfined domain walls. In the deconfined phase a single quark in 
a cubic C-periodic volume has a finite free energy. In a long cylindrical volume, 
on the other hand, the free energy increases with the length of the cylinderc3. 
Hence, the quark is "confined" even in the deconfined phase, despite the fact 
that there is no flux string emanating from the quark. "Confinement" in the 
deconfined phase is a consequence of distinct deconfined phases being aligned 
along the long direction of the cylinder. This is an indication that spontaneous 
center symmetry breaking is not just a Euclidean artifact. The sensitivity of 
the quark's free energy to the shape of the spatial volume is observable in 
Minkowski space-time and is characteristic for a spontaneously broken Z(3) c 
symmetry. Other interesting subtleties of the center symmetry are related to 
the way it arises in gauge- fixed formulations of the theoryED and to the 't Hooft 
loop order parameter c3. 



1.4 Organization of the Rest of the Article 

The rest of this article contains a detailed and more mathematical discus- 
sion of the topics outlined in the introduction. Section 2 describes the bulk 
physics of gluons at high temperature. In particular, effective theories for the 
Polyakov loop are constructed at non-zero temperature and at non-zero chem- 
ical potential for static quarks. The nature of confined and deconfined phases 
is discussed and the transitions between them are characterized. In section 3 
we study the properties and the dynamics of interfaces separating coexisting 
phases of confined and deconfined gluons and gluinos. In particular, the uni- 
versal phenomenon of complete wetting is explored using effective theories for 
the Polyakov loop. Section 4 deals with static quarks at finite volume and the 
paradoxical "confinement" in the deconfined phase. Finally, section 5 contains 
our conclusions. 
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2 Gluons at High Temperatures 



In this section we introduce SU(N) pure Yang-Mills gauge theory at non-zero 
temperature. We discuss the problem in the language of continuum field theory 
but ignore the subtleties that arise in the regularization of ultraviolet diver- 
gences. Consequently, the resulting expressions are sometimes rather formal. 
We are interested in non-perturbative questions that can in general not be 
fully addressed directly in the continuum. Using the lattice regularization, it 
is straightforward to replace the formal quantities by completely well defined 
mathematical expressions. We prefer not to do this here because mathematical 
rigor might obscure the basic physics that can be well understood in continuum 
field theory language. Still, we like to emphasize that our understanding of 
the center symmetry and its spontaneous breakdown is to a large extent due 
to studies in lattice field theory — the most powerful tool for understanding 
non-perturbative field theory from first principles. 

2.1 Path Integral and Euclidean Action 

Let us consider a system of gluons at a non-zero temperature T = 1/(3. Its 
partition function is given by 



Here H is the Yang-Mills Hamilton operator and P is a projection operator 
that enforces the Gauss law. It reduces the trace over all states to those that 
are gauge invariant and hence physical. We note that exp(— f3H) is similar 
to the time-evolution operator exp(— iHt), except that (3 plays the role of an 
imaginary (or Euclidean) time interval it. However, in order to motivate path 
integrals in Euclidean time, there is no need to start out in real Minkowski 
time and then Wick-rotate. Our starting point is equilibrium thermodynamics 
which makes no reference to real time evolution and puts us directly into Eu- 
clidean space-time. One can write a path integral expression for the partition 
function, 



Z = Tr [exp(-/3H)P]. 



(1) 




(2) 



where 





is the Euclidean action. Here g is the gauge coupling and 



Fp,v — d^A u d u A^ + [Ap, A u ] 



(4) 
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is the field strength resulting from the anti-Hcrmitian non-Abelian vector po- 
tential 

A„ - iA£A«. (5) 

The SU(N) gluon field is described by N 2 — 1 real- valued functions A^ (x) (with 
a G {1,2, ...,N 2 — 1}) of the Euclidean space-time points x = (x,X4). The A a 
are the ./V 2 — 1 Hermitian generators of the SU(N) algebra in the fundamental 
representation. As a consequence of the trace in eq.(|l|), the measure DA of 
the path integral contains fields that obey periodic boundary conditions, 

A^ (x, x 4 + /3) = A^ (x, x 4 ) , (6) 

in the Euclidean time direction. In the path integral the projection operator P 
that enforces the Gauss law manifests itself as the integration over the time- like 
vector potentials A 4 . 



2.2 Gauge Invariance and Center Symmetry 

The action of eq. (|3|) is invariant under gauge transformations 

VA,=g(A,+d,)g\ (7) 

with g(x) £ SU(N), under which the field strength transforms as 

9F^=gF^. (8) 

In order to maintain the boundary condition eq.(|6|) for the vector potential, 
we first consider gauge transformations that are strictly periodic in Euclidean 
time 

g(x,x A + f3) =g(x,x A ). (9) 

Every physical quantity must be invariant under these gauge transformations. 
However, in addition there are topologically non-trivial transformations El that 
are periodic up to a constant twist matrix h £ SU(N) 

g(x,X4 + 0) = hg{x,Xi). (10) 

When such a transformation is applied to a strictly periodic vector potential 
A^ it turns into 

9 A^(x,X4 + (3) = g{x,x 4 + f3)(A l _ l (x,x 4 + (3) + d l _ l )g{x,X4 + fi) 1 
= hg(x,x i )(A ll {x,x i ) + 9 /i )g(f,x 4 ) t /i t 
= h 9 Af,(x,x 4 )h\ (11) 
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The gauge transformed vector potential 9 A ll obeys the boundary condition 
eq.(g) only if 

h 9 A^x^ = 9 ^(f,ar 4 ), (12) 

which is the case only if h commutes with 9 Ay,. This limits us to twist matrices 
h in the center Z(iV) c of the gauge group SU(N). The elements of the center 
commute with all group elements and are multiples of the unit matrix 

h = zl, z = exp(2itin/N), n e {1, 2, N}. (13) 

It is important to note that the Z(iV) symmetry gets explicitly broken in the 
presence of dynamical fields that transform in the fundamental representation 
of SU(N). For example, quark fields transform as 

% = gV, (14) 

and are anti-periodic in Euclidean time, 

*(^,ai4 + /?) = -^>{S,Xi). (15) 

Hence, under a twisted transformation they transform into 

%(x,x 4 + /9) = g(x, x 4 + /3)ty(x. x 4 + /3) = -zg(x, x 4 )^(x, x 4 ) = -z 9 *(f,x 4 ). 

(16) 

To maintain the boundary condition eq. (j^) one must restrict oneself to z = 1 
so that the center symmetry disappears. 

The twisted transformations of eq.(|lO|) represent a global Z(JV) C symmetry 
of the action of eq.(||). Although the action is 7.(N) C invariant, other physical 
quantities need not necessarily be invariant under the center symmetry trans- 
formations. In particular, unlike the local SU (N) gauge symmetry, the global 
Z(iV) c symmetry can break spontaneously and should not be considered as 
a subgroup of the local gauge group. It should be pointed out that the cen- 
ter symmetry is a symmetry of the Euclidean action, not of the Hamiltonian. 
While the Hamiltonian acts at a given instant in time, the ~L{N) C transforma- 
tions are characterized by non-trivial boundary conditions in Euclidean time 
and thus affect the entire time evolution. Hence, it is meaningless to ask how 
states in the physical Hilbert space transform under the center symmetry. In- 
stead, the center symmetry is a symmetry of the spatial transfer matrix - 
the analog of the Hamiltonian that describes the evolution of the system in 
a spatial direction. The spatial transfer matrix also acts on a Hilbert space 
(however, not on the usual one) and the states in that space have definite trans- 
formation properties under ~L{N) C . In fact, at non-zero temperature the spatial 
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transfer matrix is at least as useful as the Hamiltonian. While the spectrum of 
the Hamiltonian is by definition independent of temperature, the spectrum of 
the transfer matrix is sensitive to the temperature and, for example, contains 
valuable information about phase transitions. 



2.3 The Polyakov Loop as an Order Parameter for Deconfinement 

As we have seen, dynamical fields that transform in the fundamental represen- 
tation break the Z(7V) C symmetry explicitly. Here we limit ourselves to cases 
where such fields are not present. Still, even then one can use static funda- 
mental charges to probe the physics of the adjoint dynamical fields. Static fun- 
damental charges (infinitely heavy test quarks) are described by the Polyakov 
loop (a Wilson loop closed around the periodic Euclidean time direction) Em 



= Tr "Pcxp( / dxi A 4 (x,Xa)). 
Jo 



(17) 



Here V denotes path ordering of the exponential. The Polyakov loop is a com- 
plex scalar field that depends on the spatial position x of the static color source. 
In the special case of two colors (N = 2) the Polyakov loop is always real. The 
Polyakov loop transforms non-trivially under twisted Z(iV) c transformations 



f 

9 $(f) = TrVexp( dx 4 9 A 4 (x, x A )) 
Jo 



Tr [g(x, X4 + (3)V exp( / dx 4 A^x, X/^))g[x, x^y] 



r/3 

= Tr [zg(x, x^fP exp( / dx± A±{x, X4))g(x, X4,y] = z$(x). (18) 
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On the other hand, it is invariant under strictly periodic gauge transformations 
(with z — 1) as it should as a gauge invariant physical quantity. 

The partition function of a system of gluons in the presence of a static 
infinitely heavy test quark is given by 

Z Q = J DA $(z)exp(-S[A]). (19) 

The thermal expectation value of the Polyakov loop, 

(*> = !/ DA $(£) exp(-5[A]) = ^ = exp(-/?F), (20) 
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is the ratio of the partition functions of the gluon systems with and without 
the external color source and hence measures the free energy F of the external 
static quark. At low temperatures color is confined and the free energy of a 
single quark is infinite (F = oo). Consequently, in the confined phase ($) = 0. 
At high temperatures, on the other hand, asymptotic freedom suggests that 
quarks and gluons become deconfined. Then F is finite and (<£>) =$0^0 hi 
the deconfined phase. Since $ transforms non-trivially under center symme- 
try transformations, a non-zero expectation value $o implies that the Z(iV) c 
symmetry is spontaneously broken in the high-temperature deconfined phase. 
It is unusual that a symmetry that is intact at low temperatures gets sponta- 
neously broken at high temperatures. However, we should keep in mind that 
the center symmetry is a symmetry of the spatial transfer matrix — not a 
symmetry of the Hamiltonian. In fact it is easy to understand why the center 
symmetry must break spontaneously in the high-temperature limit j3 — > 0. In 
this limit the integral in eq.([l7|) extends over shorter and shorter Euclidean 
time intervals and hence <&o — > Trl = N. 



2.4 Effective Actions for the Polyakov Loop 

To investigate the deconfinemcnt phase transition it is instructive to consider 
effective theories for the Polyakov loop. In principle, one can imagine to inte- 
grate out the gluon fields and derive the effective action for the Polyakov loop 
directly from the underlying Yang-Mills theory 

exp(-£[$]) = J DA 5[<P{x)-TrVex-p(J^ dx 4 A 4 (x,x 4 ))]exp(-S[A]). (21) 

Here the 5-functional ensures that <E> obeys eq.(p"7|). The resulting effective 
action is nonlocal and impossible to compute in practice. Fortunately this is 
not even necessary. To gain some qualitative insight into the nature of the 
deconfinement phase transition it is sufficient to construct simple effective the- 
ories that have the desired symmetry properties and to study their properties. 
This approach is particularly fruitful when one wants to study a universal 
phenomenon that is common to the Yang-Mills theory and to much simpler 
models. All universal aspects of this phenomenon can then be studied in a 
simple effective theory for the Polyakov loop. Universal behavior is naturally 
associated with second order phase transitions. Then a physical correlation 
length becomes infinitely large and the short-distance details of the dynamics 
become irrelevant. What matters are only the dimensionality of space and the 
symmetry properties of the order parameter. 

Let us first construct a simple effective theory for the Polyakov loop in 
SU(2) pure gauge theory. Because the trace of an SU(2) matrix is always 
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real, the Polyakov loop $(x) is then a real-valued field in three dimensions. 
It transforms non-trivially under the Z(2) c center transformations as ff $ = z$ 
with z = ±1. Consequently, it is natural to consider the action 

S[$] = J d 3 x + V($)] (22) 

with a Z(2) c invariant potential 

= a$ 2 + c$ 4 . (23) 

Stability of the problem requires c > 0. For a > the potential has a single 
minimum at $ = that corresponds to the low-temperature confined phase. 
For a < 0, on the other hand, there are two degenerate minima at <j> = ±<I>o = 
±^/a/2c corresponding to two high-temperature deconfined phases. The phase 
transition happens at a = and then $o — > indicating a second order phase 
transition. Hence, for SU(2) pure gauge theory one indeed expects universal 
behavior at the deconfinement phase transitionaQ. A simple lattice model with 
the same spatial dimension and the same Z(2) symmetry is the 3-dimensional 
Ising model. Indeed, detailed lattice gauge theory studies have shown that 
the deconfinement phase transition of the SU(2) Yang-Mills theory is second 
order El and in the universality class of the 3-d Ising model 0. Since the order 
parameter vanishes continuously, at a second order phase transition the various 
phases become indistinguishable at the transition. The lattice calculations also 
give the ratio of the critical temperature T c and ±he square root of the zero- 
temperature string tension a as T c /^/a = 0.69(1) El. 

For three colors (N = 3) the situation is different. First of all, the Polyakov 
loop <f> = $i + i<i>2 is now a complex scalar field with real part $i and imaginary 
part $2- Consequently, besides the Z(3) center symmetry, we now must also 
consider the charge conjugation symmetry!] Charge conjugation replaces the 
non-Abelian vector potential by its complex conjugate 

% = A;. (24) 

Accordingly the Polyakov loop also turns into its complex conjugate 

c $(f) = TrVexp( dx A c A 4 (x, x A )) 
Jo 

r& 

= TrPexp(/ dx± A 4 (x,x 4 )*) = $(£)*. (25) 



''For two colors charge conjugation is equivalent to a global gauge transformation. 
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The natural candidate for an effective action now is 

5[*] = J d 3 x [^di&di® + Vm , (26) 

with a Z(3) c and charge conjugation invariant potential, i.e. 

V(z<f>) = V(<Z>*) = V($). (27) 

The most general quartic potential that obeys these symmetries is given by 

V($) = a|$| 2 + 6$i(3$2 - $?) + c|$| 4 . (28) 

Again, in order to have a stable system we need c > 0. For a > the system 
has a minimum at $i = $2 = corresponding to the confined phase. For 
< a < 96 2 /32c there is another minimum at 

, , 3b + V% 2 ~32ac 

$1 = $ = , $ 2 = 0, 29 

8c 

which has two other Z(3) c copies. We denote the four minima by 
$M = $ , $( 2 ) =_(_! + iV3)$ , * (3) = ^(-1 - «V3)$o, $ (4) = 0. (30) 

The first three minima correspond to the three high-temperature deconfined 
phases. For b 2 = Aac all four minima are degenerate. This corresponds to the 
phase transition temperature. However, in contrast to the SU(2) case, at the 
phase transition we now have $0 = \/a/c = b/2c ^ which corresponds to 
a first order phase transition. At a first order phase transition one does not 
have universal behavior in the bulk because the bulk correlation length remains 
finite and the short distance details of the dynamics remain important at the 
phase transition. A simple lattice model in the same spatial dimension and 
with the same global symmetries as the SU(3) pure gauge theory is thepi- 
d 3-state Potts model u which indeed has a first order phase transition 113. 
Furthermore, detailed numerical studies in lattice gauge theory have shown 
that the deconfinement phase transition in SU(3) Yang-Mills theory is also 
first orderll 2 ]. In this case the ratio of the critical temperature and the square 
root of the string tension is T c / 'y/a — 0.64(1). 

2. 5 Gluons in a Background of Static Quarks 

As we have seen, the partition function of a system of gluons in the presence 
of a single static quark at position x is given by eq. (|l9|) . Similarly, the gluon 
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partition function in the presence of a single anti-quark is 



Zq = j DA $(£)*exp(-S[4]). (31) 

Consequently, the partition function of a system of gluons in the background 
of n static quarks and n static anti-quark takes the form 



z n .n = J DA ^[J d " x ^ n k [ J d3y Hy)TeMS[A]). 



(32) 



Here we have integrated over the positions x of quarks and y of anti-quarks. 
The factors n\ and n\ take into account that quarks as well as anti-quarks are 
indistinguishable particles. We consider quarks of mass M with quark chemical 
potential \i that couples to the baryon number B = (n — n)/3 and thus obtain 

Z(ji) = ^Z njr exp(-/?M(n + n))cxp(/? M (n-n)) 

= J DAj2^{j d3x $ (^)]" exp(-/3(M - M )n) 

x ^ifd 3 y $(^)Texp(-/3(M + M )n)exp(-S[A]) 

= J L>^cxp[ J d 3 x exp(-/3(M - /i))] 

x cxp[ J d 3 y <t>(y)* exp(-/3(M + //))] exp(-S[A]). (33) 

Consequently, the effective action for the gluons now takes the form 

S eff [A] = S[A]-J d 3 x [$(f)exp(-/3(M- M )) + $(f)*exp(-/3(M+/i))]. (34) 

As expected in the presence of fundamental charges, the additional terms in 
the effective action break the center symmetry explicitly. In the presence of 
the chemical potential, the effective action is in general complex. The SU(2) 
case is exceptional because then the Polyakov loop and hence the effective 
action take real values. The action is real even in the SU (N) case for purely 
imaginary values of the chemical potential. In that case, /i = iB^ enters the 
action as a constant vector potential B4 in the 4-direction that couples to the 
baryon number. A real chemical potential, on the other hand, enters the action 
as a purely imaginary vector potential. Taking the limit M, jj, — > 00 with M — fi 
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finite we obtain 



S ef f[A] = S[A] - J d 3 x exp(-/3(M - /*)). (35) 

The chemical potential favors positive real values of 

In the simple $ 4 model for the SU(3) Polyakov loop the chemical potential 
enters the effective potential as 

V($) = a|$| 2 + 6$i(3$2 - $1) + c|$| 4 - W&, (36) 

where ft. represents exp(— /?(M — /i)). At non-zero h the deconfined phase d^ 
with a real valued Polyakov loop is favored while the other two deconfined 
phases are suppressed. The deconfined phase d^' coexists with the confined 
phase along a first order transition line that terminates in a critical endpoint. 
Since the Z(3) c symmetry is explicitly broken in the presence of the chemical 
potential, at non-zero h the Polyakov loop is non-zero even in the confined 
phase. We denote the values of the Polyakov loop in the confined and the 
deconfined phase by and <&( d \ respectively. Along the first order transition 
line at which the confined and the deconfined phase coexist, the potential (at 
$2 = 0) can then be written as 

V(&) = c(*i - $W) a ($i - $ (c) ) 2 + V . (37) 

Comparing eq.(|36j) and eq.(|37|) one derives a condition for the location of the 
first order phase transition line 

4ac 8c 2 , ,„„, 

— -l = —h. (38) 

It is non-trivial to analyze the phase diagram with lattice field theory methods 
because the_complex action prevents the use of standard importance sampling 
techniques cB. Using a cluster algorithm, the complex action problem has been 
solved for the 3-d 3-state Potts modeled. The corresponding phase diagram is 
shown in figure 1. In the Potts model the parameter k — k c plays the role of 
1 — Aac/b 2 . Interestingly, the transition line in the Potts model is consistent 
with a straight line as predicted by eq.(p8|). The discontinuity of the Polyakov 
loop across the transition takes the form 



It vanishes at the critical endpoint h — He in which the first order line ter- 
minates at a second order phase transition. A critical exponent (3 (not to be 
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Figure 1: The phase diagram of the 3-d 3-state Potts model in the (h, n)-plane. Qualitatively, 
k plays the role of the temperature in Yang-Mills theory. The ordinary deconfinement phase 
transition at h = is a triple point T from which a line of first order phase transitions 
emerges. This line terminates in the critical endpoint E. 

confused with the inverse temperature) is defined by 

_ $(c) K (! _ *L)£ (4 ) 
He 

The simple mean field theory result of eq.(^) yields (3 — 1/2. However, this 
result is not reliable. Instead, one expects the endpoint to be in the universality 
class of the 3-d Ising model which has f3 = 0.3260(8). A similar endpoint has 
been studied numerically by Karsch and Stickan in the 3-d 3-state Potts model 
with a real action Ell. Indeed they find very good agreement with universal 3-d 
Ising behavior. 

3 Interfaces Separating Confined and Deconfined Phases 

At a first order phase transition the various phases remain distinguishable 
because the order parameter is discontinuous at the transition ($o 7^ 0). Con- 
sequently, at a first order phase transition distinct phases coexist with one 
another. The phases are spatially separated from each other by interfaces with 
a non-zero free energy per unit area (the interface tension) . In SU (3) Yang- 
Mills theory there are confined-deconfined interfaces that separate the low and 
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high-temperature phases which coexist at the first order deconfinement phase 
transition. In SU(2) gauge theory, on the other hand, such interfaces do not 
exist because the confined and the deconfined phase become indistinguishable 
at the second order phase transition. Both in SU(2) and SU(3) Yang-Mills 
theory there are deconfined-deconfined domain walls that separate the vari- 
ous deconfined phases at high temperatures. When a deconfined-deconfined 
domain wall in the SU(3) theory is cooled down to the phase transition, it 
splits into a pair of confined-deconfined interfaces. Then the confined phase 
forms a complete wetting layer that grows to a macroscopic size at the phase 
transition. At the phase transition the size of the wetting layer diverges and 
again the microscopic details of the dynamics become irrelevant. 

3.1 Deconfined- Deconfined Domain Walls and Confined- Deconfined Interfaces 

When a discrete symmetry gets spontaneously broken, there exist a number of 
degenerate phases related to one another by symmetry transformations. These 
phases can coexist and are separated by domain walls. In the gluon plasma 
the spontaneous breakdown of the Z(3) c center symmetry gives rise to three 
distinct deconfined phases. These phases are separated from one another by 
deconfined-deconfined domain walls. The domain walls are stable at all tem- 
peratures above T c . Due to asymptotic freedom, the domain wall tension (free 
energy per unit area) can be computed semi-classically in the high-temperature 
limit. The domain wall is a soliton in the effective potential for the Polyakov 
loop. The 1-loop effective potential for the complex phase of the Polyakov loop 
was first calculated by Weiss E3. This calculation shows that there are three 
degenerate minima corresponding to three distinct deconfined phases. Based 
on this result, Bhattacharya, Gocksch, Korthals Altes and Pisarski have ana- 
lytically determined the_tension of a deconfined-deconfined domain wall in the 
high-temperature limit Ej 

a dd = ^T 3 , (41) 

where g is the gauge coupling renormalized at the scale T. 

Interfaces separating coexisting phases also occur at any first order phase 
transition. In that case, the coexisting phases are in general not related by 
symmetry. For example, in SU(3) Yang-Mills theory the high-temperature 
gluon plasma coexists with the low-temperature glueball phase at T c . The 
corresponding confined-deconfined interfaces are stable only at T c . Unlike the 
deconfined-deconfined domain wall tension, the confined-deconfined interface 
tension cannot be computed semi-classically since T c is too small for pertur- 
bation theory to be applicable. Instead there have been several attempts to 
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extract the interface tension from lattice calculations It should be noted 
that these calculations do not contain dynamical quark effects. Due to rather 
strong lattice spacing artefacts, it is difficult to reliably extrapolate the results 
to the continuum limit. A rough estimate of the confined-deconfined interface 
tension gives a c d ~ 0. 020(5)7^. This result is interesting from a cosmological 
point of view. The confined-deconfined interface tension sets a scale for the 
nucleation rate of bubbles of confined phase from the quark-gluon plasma and 
hence for the spatial inhomogeneities generated at temperatures around 100 
MeV. These inhomogeneities may have an effect on the primordial synthesis 
of light atomic nuclei which takes place when the universe has cooled down 
to temperatures in the 1 MeV range. The small value for a c d from above 
suggests that the inhomogeneities are small and that the standard scenario of 
primordial nucleosynthesis is correct. 



3.2 Complete Wetting of Hot Gluons 

Frei and Patkos were first to conjecture that complete wetting occurs in ST/ (3) 
pure gauge theory at the high-temperature deconfinement phase transition 
EH. The interface tensions a c d and add of confined-deconfined interfaces and 
deconfined-deconfined domain walls, respectively, determine the shape of a 
droplet of confined phase that wets a deconfined-deconfined domain wall. As 
shown in figure 2a, such a droplet forms a lens with opening angle 9, where 

9 . N 

a d d = 2a cd cos -. (42) 

This equation for 9 follows from the forces at the corner of the lens being in 
equilibrium. In condensed matter physics the inequality 

add < 2a cd , (43) 

was derived by WidomEl. It follows from thermodynamic stability because a 
hypothetical deconfined-deconfined domain wall with add > 2a c( j would simply 
split into two confined-deconfined interfaces. If add < 2a C£ ; the confined phase 
forms droplets at a deconfined-deconfined domain wall. In condensed matter 
physics this phenomenon is known as incomplete wetting. For add = 2ct c d, on 
the other hand, 9 — and the lens shaped droplet degenerates to an infinite 
film, as shown in figure 2b. When such a film is formed, this is called complete 
wetting. 

Complete wetting is a universal phenomenon of interfaces characterized by 
several critical exponents. For example, the width r of the complete wetting 
layer diverges as 

roc(T-T c )-*, (44) 
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Figure 2: Incomplete versus complete wetting, (a) For ct^d < 2ct c( j one has incomplete 
wetting with 9^0. Then the confined phase forms a lens shaped droplet at the deconfined- 
deconfined domain wall, (b) Complete wetting corresponds to a^d = 2« c( j- Then 9 = 0, and 
the confined phase forms a film that splits the deconfined- deconfined domain wall into two 

confined- deconfined interfaces. 



where T — T c measures the deviation from the point of phase coexistence and ip 
is a critical exponent. The value of ip depends on the range of the interactions 
between the two interfaces that enclose the complete wetting layer. In SU(3) 
Yang-Mills theory the interactions between the interfaces are short-ranged and 
the interaction energy per unit area is given by 7exp(— r/ro). In addition, 
above T c the confined phase that forms the wetting layer has a slightly larger 
bulk free energy than the deconfined phases. Close to T c , the additional bulk 
free energy per unit area for a wetting layer of width r is given by 6(T — T c )r. 
Hence, the total free energy per unit area of the two interface system relative 
to the free energy of two infinitely separated interfaces at T c is given by 

a dd (T) - 2a cd {T c ) = 7 exp(-r/r ) + 5(T - T c )r. (45) 

Minimizing the free energy with respect to r, one finds the equilibrium width 

r = -ro log ^(T-T c ), (46) 
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such that ijj = 0. Inserting the equilibrium value of r one finds 

a dd (T)-2a cd (T c ) = 6r (T-T c )[l-log^(T-T c )]. (47) 

7 

This-is consistent with the critical behavior in the Z(3) c symmetric 3> 4 the- 
orylla. As we will see later, the same critical exponents follow for the high- 
temperature deconfinement transition in the supersymmetric case. 

The universal aspects of the dynamics of confined-deconfincd interfaces 
and deconfined-deconfined domain walls in SU(3) Yang-Mills theory_can be 
investigated using the effective action for the Polyakov loop of cq. (pq) E3. The 
corresponding potential takes the form V(&) = a|$| 2 + M>i(3$2 — $j) + c|<I>| 4 of 
eq.(|2^). The deconfinement phase transition occurs at b 2 — Aac because then 
the confined phase minimum at $*- 4 ^ =0 is degenerate with the deconfined 
phase minimum at = <J> an d its Z(3) c copies at $( 2 ) and $( 3 ). The profile 
$(£) = &(z) of a planar interface perpendicular to the z-direction follows from 
the solution of the corresponding classical equations of motion 

2a$i + 3b(<f>l - $?) + 4c$i($ 2 + $§), 

2a$ 2 + 66$i$2 + 4c$ 2 (*i + $2)- ( 48 ) 

The solution for a confined-deconfincd interface interpolating between the con- 
fined phase at negative z and the deconfined phase <I>( 1 - ) at positive z takes the 
form 

$ x (z) = i*„ [1 + tanh^z)], $ 2 (z) = 0. (49) 

Here [i — bj \/8c determines the width of the interface. The confined-deconfincd 
interface tension a cd is given by 

f°° 1 f°° 
Pa cd = / dz [-d&di* + V($)] =2 dz V($) = f^. (50) 

J — oo * J — oo 

The interface profile, its width, as well as its tension are non-universal prop- 
erties that are quantitatively different in the simple effective <& 4 theory and in 
SU{3) Yang-Mills theory. 

Still, there are universal properties related to the phenomenon of complete 
wetting. To investigate this, we now consider the profile of a deconfined- 
deconfined domain wall that interpolates between the deconfined phases 
and $' 3 ^. The domain wall profile is a solution of eq.([48|) with the condition 
that $ approaches <I>( 2 ) at z = oo and $' 3 ' at z = — oo. A typical solution 
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Figure 3: Profile of a deconfined-deconfined domain wall in SU(3) Yang-Mills theory deep 

in the deconfined phase. 



deep in the deconfined phase is illustrated in fig. 3. There is a sharp transition 
from one deconfined phase to the other. Figure 4 shows the domain wall profile 
very close to the deconfinement phase transition. In this temperature region 
the deconfined-deconfined domain wall splits into a pair of confincd-deconfined 
interfaces that enclose a complete wetting layer of confined phase. Very close 
to the phase transition one can even construct an analytic solution of the 
completely wet deconfined-deconfined domain wall by combining two solutions 
for confined-deconfincd interfaces 

&i(z) = — i$ [2 + tanh^i(z — z ) — tanh/i(z + z )], 
V3 

$ 2 (z) = — $ [tanh^(z-z )+tanh y u(z + z )]. (51) 

The interface separation and hence the width of the confined complete wetting 
layer is given by 

1 , .1 2ac, 
^2z = --log(---). (52) 

Since the phase transition corresponds to b 2 = Aac, the width of the wetting 
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Figure 4: Profile of a deconfined-deconfined domain wall in SU(3) Yang-Mills theory close 
to the phase transition where the wall splits into two confined- deconfined interfaces with a 
complete wetting layer of confined phase between them. 



layer diverges logarithmically at T c . This is consistent with the critical ex- 
ponent if) = as expected for interfaces with short-range interactions. The 
interface tension of the deconfined-deconfined domain wall is given by 

^ = ^-2(1-^)]. (53) 

Hence, consistent with the leading (T — T c ) dependence in eq.(|47|), one obtains 

p(2a cd - a dd ) = 2^1(1 - -^). (54) 

3.3 Strings Ending on Walls in Super 'symmetric Yang- Mills Theory 

Witten has argued that color flux strings in supersymmetric Yang-Mills theory 
can end on domain walls El This effect follows from a calculation in M-theory, 
where the domain wall is represented by a D-brane on which strings can end. 
In supersymmetric Yang-Mills theory a Z(3) x chiral symmetry an unbroken 
remnant of the anomalous U(1)r symmetry — is spontaneously broken at low 
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temperatures by a non-zero value of the gluino condensate \- As a consequence, 
there are three distinct confined phases, characterized by three different values 
of x which are related by Z(3) x transformations. Regions of space filled with 
different confined phases are separated by confined-confined domain walls. The 
properties of domain walls in supersymmetrinjthcories have been investigated in 
great detail by_Shifman and collaboratorsHo and also by Carroll, Hellerman 
and Trodden c3. Often, topological defects which result from spontaneous 
symmetry breaking have the phase of unbroken symmetry at their cores. For 
example, magnetic monopoles or cosmic strings have symmetric vacuum at 
their centers. Here we construct an effective theory for the Polyakov loop 
and the gluino condensate in which Z(3) x restoration implies the breaking 
of the Z(3) c center symmetry. Then the high-temperature deconfmed phase 
appears at the center of the domain wall forming a complete wetting layer. 
As a consequence, the Polyakov loop has a non-zero expectation value there 
and thus a static quark has a finite free energy close to the wall so that its 
string can end there. In supersymmetric Yang-Mills theory strings can end 
on confined-confined domain walls because they can transport color flux to 
infinity at a finite free energy cost. A similar scenario iruiolving a non-Abelian 
Coulomb phase has been discussed at zero temperature L.1 

In supersymmetric SU(3) Yang-Mills theory the Z(3) x chiral symmetry 
is spontaneously broken in the confined phase. The corresponding order pa- 
rameter is the complex valued gluino condensate \ — Xi + *X2- Under chiral 
transformations z G Z(3) x the gluino condensate transforms into x' — X z 
and under charge conjugation it gets replaced by its complex conjugate. At 
high temperatures one expects chiral symmetry to be restored and — as in the 
non-supersymmetric theory — the Z(3) c center symmetry to be spontaneously 
broken due to deconfinement. Consequently, the effective action describing the 
interface dynamics now depends on both order parameters $ and x an( i 

S[*,X\ = / d 3 x [^a i $+i%*9 iX +F($,x)]. (55) 

The most general quartic potential consistent with Z(3) c , Z(3) x and charge 
conjugation symmetry now takes the form 

V($,X) = a|$| 2 + 6$i(3<^-$?) + c|<I>| 4 

+ rf|x| 2 + exi(3x^Xi) + /lx| 4 + .9l'i>| 2 |x| 2 . (56) 

Assuming that deconfinement and chiral symmetry restoration occur at the 
same temperature and that the phase transition is first order, three chirally bro- 
ken confined phases coexist with three distinct chirally symmetric deconfmed 
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Figure 5: Profile of a confined- confined domain wall in supersymmetric SU(3) Yang-Mills 
theory deep in the confined phase. Note that $1(0) ^ 0, i.e. the center of the wall has 
properties of the deconfined phase. 



phases. The three deconfined phases have = <J> , = (-l/2+iv / 3/2)$o 
and $( 3 > = (-1/2-^/2)^0 and x (1) = X (2) = X (3) = 0, while the three con- 
fined phases are characterized by $^ 4 -* = <&( 5 - ) = < J> < - 6 ) = and x'' 4 '' = Xo <= R> 
X (5) = (-l/2 + iV3/2)xo and x ( ' 6) = (-1/2 - n/3/2)xo- The phase transition 
temperature corresponds to a choice of parameters a,b,...,g such that all six 
phases <£>("), ^(") represent degenerate absolute minima of V($,x). 

We now look for solutions of the classical equations of motion, representing 
planar domain walls, i.e. $(x) = $(2), x(^) = x( z )i where z again is the 
coordinate perpendicular to the wall. The equations of motion then take the 
form 

f^^m^ fx 1= dv 

dz 2 dfy' dz 2 d Xl ' 

Figure 5 shows a numerical solution of these equations for a domain wall 
separating two confined phases of type and i.e. with boundary con- 
ditions $(00) = <f>( 5) , x(°°) = X^ and $(-00) = x( _ °°) = x^ ■ Figure 
5 corresponds to a temperature deep in the confined phase. Still, at the do- 
main wall the Polyakov loop is non-zero, i.e. the center of the domain wall 
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Figure 6: Profile of a confined- confined domain wall in supersymmetric SU(3) Yang-Mills 
theory close to the phase transition. In this region the wall splits into two confined- deconfined 
interfaces with a complete wetting layer of deconfined phase between them. 



shows characteristic features of the deconhncd phase. Figure 6 corresponds to 
a temperature very close to the phase transition. Then the confined-confined 
domain wall splits into two confined-deconfined interfaces and the deconfined 
phase forms a complete wetting layer between them. The solutions of figure 5 
have deconfined phase of type S 1 ^ at their centers. Due to the Z(3) c symmetry, 
there are related solutions with deconfined phase of types S 2 ^ and . 

For the special values d = a = 0, e = b, f = c, g = 2c one can find an 
analytic solution for a confined-deconfined interface. Combining two of these 
solutions to a confined-confined domain wall, one obtains 

$i0) = -i$ [tanh/i(z-z ) -tanh/i(z + z )], $ 2 (z) =0, 

Xi(z) = ~^Xo[2 + tanh^O - z ) - tanh^(z + z )], 

X2O) = — xoNnh^O - z ) + tanh^fz + z )}, (58) 

where $0 = Xo = 3b/ 4c and jj, = 3b/4^fc. The critical temperature corresponds 
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to e 4 // 3 = b 4 /c 3 . Near criticality, where A = e 4 // 3 — 6 4 /c 3 is small, the above 
solution is valid up to order A 1 / 2 , while now e = b and / = c are satisfied to 
order A. The width of the deconfined complete wetting layer, 

r = 2z = -^logA + C, (59) 

where C is a constant, grows logarithmically as we approach the phase tran- 
sition temperature. This is the expected critical behavior for interfaces with 
short-range interactions, which again implies the critical exponent ip = 0. 

Now we wish to explain why the appearance of the deconfined phase at 
the center of the domain wall allows a QCD string to end there. We recall that 
an expectation value ($) ^ implies that the free energy of a static quark is 
finite. Indeed, the solution of figure 5 containing deconfined phase of type 
has $i(0) so that a static quark located at the center of the wall has finite 
free energy. As one moves away from the wall, the Polyakov loop decreases as 

$i(z) cx cxp(-/3F(z)) cx cxp(-z v /2a + 2g X 2 ). (60) 

Consequently, as the static quark is displaced from the wall, its free energy 
F{z) increases linearly with the distance z from the center, i.e. the quark is 
confined to the wall. The string emanating from the static quark ends on the 
wall and has a tension 

a=lim E^l = lj2a + 2g X l (61) 

z^oo z p v 

Still, there are the other domain wall solutions (related to the one from above 
by Z(3) c transformations), which contain deconfined phase of types S 2 ^ and 
One could argue that, after path integration over all domain wall con- 
figurations, one obtains ($) = 0. However, this is not true. In fact, the 
wetting layer at the center of a confined-confined domain wall is described by 
a two-dimensional field theory with a spontaneously broken Z(3) c symmetry. 
Consequently, deconfined phase of one definite type spontaneously appears at 
the domain wall. 

4 Paradoxical "Confinement" in the Deconfined Phase 

We now want to study some subtle effects related to the spontaneous break- 
down of the center symmetry at high temperatures. Whenever a global sym- 
metry gets spontaneously broken, the system becomes very sensitive to the 
spatial boundary conditions. To expose this sensitivity explicitly, it is often 
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very instructive to work in a finite spatial volume. This is much better than 
working in an infinite volume with unspecified boundary conditions, which 
tends to obscure the physics in the infrared. Instead, by consistently working 
in a finite volume, we learn how the infinite volume limit can be approached 
in a meaningful way. 

4-1 Periodic Versus C-Periodic Spatial Boundary Conditions 

Most numerical studies of lattice gauge theories use periodic spatial boundary 
conditions, 

Ap(x + Ltei,Xi) = An(x,X4,). (62) 

Here Li is the size of a 3-dimensional torus in the i-direction and ej is the 
unit-vector pointing in that direction. However, periodic boundary conditions 
are not the most useful ones for studying a single Polyakov loop. As we have 
seen, a Polyakov loop represents a single static quark that transforms non- 
trivially under the center of the gauge group. Indeed, the test quark carries 
a non-zero ~L{N) C charge. For topological reasons, a single charged particle 
cannot exist in a periodic volume. The flux emanating from the charge cannot 
escape to infinity and conseauently the Gauss law demands that it must end in 
a compensating anti-chargeE3. As a whole, a periodic system is always neutral. 

Here we consider periodic boundary conditions in .thje x and y-directions 
and C-periodic boundary conditions in the z-directionE3Ell. When a C-periodic 
field is shifted by L z , it is replaced by its charge conjugate, i.e. for C-periodic 
gluons 

A^(x + L z e z ,t) = c A tl {x,t)^A f ,{x,ty. (63) 

In a C-periodic volume, a single quark can exist, because now its center elec- 
tric flux can escape to its charge conjugate partner on the other side of the 
boundary. Note that the system is still translationally invariant. The allowed 
gauge transformations, as well as the Polyakov loop, also satisfy C-periodicity 

g{x + L z e z ,t)=g{x,t)*, $(£ + L z e z ) = $(£)*. (64) 

As a consequence of the boundary conditions, the Z(3) c center symmetry is 
now explicitly broken t3. If one again considers a transformation g(x,t + (3) — 
zg(x,t), which is periodic, up to a center element z E Z(3) c , in the Euclidean 
time direction, one finds 

g(x,t)* = g{S,t + 0)*z = g(x + L z e z ,t + (3)z = g(x + L z e z ,t)z 2 =g(x,t)*z 2 . 

(65) 

Consistency requires z 2 = 1 and hence 2 = 1 (because z £ Z(3) c ). Of course, 
in the infinite volume limit, the explicit Z(3) c symmetry breaking due to the 
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spatial boundary conditions disappears. With C-periodic boundary conditions, 
($) is always non-zero in a finite volume. In the confined phase ($) goes to 
zero in the infinite volume limit, while it remains finite in the high-temperature 
deconfined phase. C-periodic boundary conditions are well-suited for studying 
the free energy of single quarks, while with periodic boundary conditions a 
single quark cannot even exist. 



4-2 Static Quarks in C-periodic Cylinders below and above T c 



First, we consider the system in the confined phase at temperatures T < T c 
in a cylindrical volume with cross section A = L x L y and length L z 3> L x , L y . 
Diagrammatically, the partition function is 



Z : 



exp{-Pf c AL 2 



(66) 



where f c is the temperature-dependent free energy density in the confined 
phase. The expectation value of the Polyakov loop (times Z), on the other 
hand, is given by 



zm = 



exp(-/3/ c AL z )E exp(-f3aL z 



(67) 



where a is the string tension, which is again temperature-dependent. The 
confining string (denoted by the additional line in the diagram) connects the 
static quark with its anti-quark partner on the other side of the C-periodic 
boundary. Hence, the free energy of the quark is given by 



F = log S 



aL, 



(68) 



The free energy diverges as L z — > oo, indicating that the quark is confined. 

Now let us consider the deconfined phase at temperatures T > T c , where 
three distinct deconfined phases coexist. In a cylindrical volume, a typical 
configuration consists of several bulk phases, aligned along the ^-direction, 
separated by deconfined-deconfined domain walls. The domain walls cost free 
energy F proportional to their area A, such that their tension is given by 
&dd = F I A. What matters in the following is the cylindrical shape, not the 
magnitude of the volume. In fact, the cylinders can be of macroscopic size. 
The expectation value of the Polyakov loop-in a cylindrical volume can be 
calculated from a dilute gas of domain walls E3. The domain wall expansion of 
the partition function can be viewed as 



Z = 



+ 






+ 







(69) 
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The first term has no domain walls and thus the whole cylinder is filled with 
deconfincd phase d^ only. An entire volume filled with either phase S 2 ^ or d^ 
would not satisfy the boundary conditions. The second and third terms have 
one domain wall separating phases d^ and d^ . Here, C-periodic boundary 
conditions exclude phase d^\ The sum of the diagrammatic terms is given by 

Z = cxp(-f3f d AL z ) 

+ 2/ dz exp(-f3f d Az)jexp(-(3a dd A)exp{-f3f d A(L z - z)) + ... 
Jo 

= exp{-(3f d AL z )[l + 2'yexp{-0a dd A)L z + ...}. (70) 

The first term is the Boltzmann weight of deconfincd phase of volume AL Z 
with a free energy density f d . The second term contains two of these bulk 
Boltzmann factors separated by a domain wall contribution 7exp(— (3a dd A), 
where 7 is a factor resulting from capillary wave fluctuations of the domain 
wall. In the above expression, we have integrated over all possible locations z 
of the domain wall. It is straightforward to sum the domain wall expansion to 
all orders, giving 

Z = exp(-/3f d AL z + 2 1 exp(-pa dd A)L z ). (71) 
In exactly the same way one obtains 

($) = $ exp(-3 7 exp(-(3a dd A)L z ). (72) 
The free energy of a static quark in a C-periodic cylinder is therefore given by 

F = -^log^ + ^-cxp(-pa dd A)L z . (73) 

This result is counter intuitive. Although we are in the deconfined phase, the 
quark's free energy diverges in the limit L z — > 00, as long as the cross section 
A of the cylinder remains fixed. This is the behavior one typically associates 
with confinement. In fact, 

a' = ^exp(-[3a dd A) (74) 

plays the role of a "string tension" , even though there is no physical string 
that connects the quark with its anti-quark partner on the other side of the 
C-periodic boundary. "Confinement" in C-periodic cylinders arises because 
disorder due to many differently oriented deconfined phases destroys the cor- 
relations of center electric flux between quark and anti-quark. 
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Of course, this paradoxical confinement mechanism is due to the cylin- 
drical geometry and the specific boundary conditions. Had we chosen C- 
periodic boundary conditions in all directions, the deconfined phases S- 2 ^ and 

would be exponentially suppressed, so that the entire volume would be 
filled with phase only. In that case, the free energy of a static quark is 
F = — log$o, which does not diverge in the infinite volume limit, as 
expected. Had we worked in a cubic volume, a typical configuration would 
have no domain walls, the whole volume would be filled with deconfined phase 
d^ and again the free energy of a static quark would be F = — log<I>o. 
Finally, note that, in a cylindrical volume, even the static Coulomb potential 
is linearly rising with a "string tension" e 2 /A. Due to Debye screening, this 
trivial confinement effect is absent in the gluon plasma. 

The observed paradoxical confinement implies that deconfined-deconfincd 
domain walls are more than just Euclidean field configurations. In fact, they 
can lead to a divergence of the free energy of a static quark and thus they 
have physically observable consequences even in Minkowski space-time. Of 
course, the issue is rather academic. First of all, the existence of three dis- 
tinct deconfined phases relies on the Z(3) c symmetry, which only exists in a 
pure gluon system — not in the real world with light dynamical quarks. Sec- 
ond, the effect is due to the cylindrical geometry and our choice of boundary 
conditions. Even though it is not very realistic, this set-up describes a per- 
fectly well-defined Gedanken experiment, demonstrating the physical reality 
of deconfined-deconfined domain walls. 

5 Conclusions 

Thanks to the efforts of numerous people, the role of the center symmetry 
for the confinement and deconfinement of static quarks is by now very well 
understood and is more or less a closed subject. In this article we have con- 
centrated on those aspects of the problem that can be easily explained using 
simple analytic calculations in continuum field theory. However, we again want 
to emphasize the very important role that lattice field theory, in particular, 
numerical simulation of SU(2) and 51/ (3) Yang-Mills theory has played in 
understanding the dynamical role of the center symmetry It is fair to say 
that lattice gauge theory has led to a reliable high-precision numerical solu- 
tion of Yang-Mills theories. Due to severe numerical problems in dealing with 
fcrmionic fields, this can presently not be said about full QCD including dy- 
namical quarks. Still, continuous progress is being made in applying lattice 
methods and it is to be expected that an accurate numerical solution of QCD 
will eventually be obtained. In contrast to Yang-Mills theory, the center sym- 
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metry plays a minor role in full QCD because it is explicitly broken in the 
presence of quarks. However, along the way, other theories for example, 
supersymmetric Yang-Mills theories — can be solved with lattice field theory 
methods as well. At that point the center symmetry may become an inter- 
esting research topic again. In any case, we hope that we have convinced the 
reader that — thanks to the center symmetry — a world without quarks would 
not be quite as boring as one might have naively expected. 
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